IFUM - 463/FT 
March 1994 



ID-Disordered Conductor with Loops 
Immersed in a Magnetic Field 



Q\ ! Alberto S. Cattaneo [] 

■ Dipartimento di Fisica, Universita degli Studi di Milano 

and INFN, sezione di Milano, 20133 Milano, Italy. 



Andrea Gamba 



Dipartimento di Matematica, Politecnico di Torino, 10129 Torino, Italy 

^ ■ and INFN, sezione di Pavia, 27100 Pavia, Italy 
O 

g ; Igor V. Kolokolov [] 

I Sezione INFN di Milano, 20133 Milano, Italy 

g ! Maurizio Martellini [] 

i-^j ■ Dipartimento di Fisica, Universita degli Studi di Milano, 20133 Milano, Italy 

Ch ! and Sezione INFN di Pavia, 27100 Pavia, Italy. 

O 
> 

X 

.s. 

Abstract 

We investigate the conductance of a 1-D disordered conducting loop with two contacts, 
immersed in a magnetic flux. We show the appearance in this model of the AFtshuler- 
Aronov-Spivak behaviour. We also investigate the case of a chain of loops distributed 
with finite density: in this case we show that the interference effects due to the presence 
of the loops can lead to the derealization of the wave function. 
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1 Introduction 



The Aharonov-Bohm effect arises combining a non-trivial topology with the presence of 
the magnetic field [Jj. When disorder is added to this picture it has been shown that 
usually observed quantities acquire new properties. For instance, the conductivity of a thin 
cylindrical conductor immersed in a magnetic flux <p is a periodic function having period <pQ, 
the quantum of magnetic flux, but in the presence of disorder one gets a period 0o/2 (see 
Ref. 0). This effect has been observed in experiments ||. In Ref. computations were 
carried on for the 2D-case in the framework of weak localization theory [Q and essential use 
was made of perturbation theory. 

The papers || [7] have arisen interest toward ID-models presenting the Al'tshuler- 
Aronov-Spivak behaviour. In these works the study of the conductivity of a ring immersed 
in a magnetic flux and having two contacts on opposite sides was carried on. The case of 
the ordered ring was studied there thoroughly. In || it was considered a disordered ring, 
but localization effects were completely ignored. 

In the present paper we consider a loop of one-dimensional disordered conductor im- 
mersed in a magnetic flux <fr (measured in units of 4>o). The geometry of the contact that we 
are considering differs from the one analyzed in Refs. ||-|§ (see section 2 and Fig. 1). This 
allows us to describe the contact in terms of glueing matrices for solutions of the Schrodinger 
equation rather than in terms of the S-matrices used in Refs. This point may turn 

useful in generalizations of the model. 

We compute the conductivity of the sample in the framework of the Landauer-Biittiker 
approach ||, |j and show that it is a periodic function of 20, so that the frequency of 
oscillations is doubled. This result is in agreement with the considerations put forth in 



Ref. M (see also pi for a review). 



2 Conductance of a loop immersed in a magnetic flux 

In order to compute the conductance of the sample depicted in Fig. 1 we resort to Landauer's 
formula ||, ||: 

where r is the transmission probability (calculated at the Fermi level), r can be obtained 
from the 2x2 transfer matrix T which maps the space of solutions ip{x) of the Schrodinger 
equation from one space point to another: 

»<*») = *■«<*.). «M=($3-£$i). '=« P) 

here k denotes the one- dimensional momentum of the electron and the Schrodinger equation 
has the form 

- ifj" + U(x)ip = k 2 i). (3) 
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Figure 1: The form of the conducting loop. 

Let us determine the transfer matrix T for a generic loop. The one-dimensional description is 
adequate only asymptotically, far away from the cross-junction domain. Inside this domain 
we should use the three-dimensional Schrodinger equation. However, if the domain is small 
enough (i.e., of size ~ 1/k), we just need to get some linear relations among the effective 
asymptotic one-dimensional wave functions. Since the Schrodinger equation is a second 
order PDE, the fixing of the boundary values of the wave function determines the solution 
uniquely with all its derivatives. 




d a 
Figure 2: The cross-junction. 

This means that among the eight quantities ip(a), ip'(a), ip(b), ip'(b), ip(c), ip r (c), ip(d), ip'(d) 
(see Fig. 2) there should exist four linear relations, which can be written for a symmetric 
junction as 

(:H* *)(::> 

with A and B being 2x2 complex matrices. If we now connect b and c through a generic 
transfer matrix To (or e^To in the presence of a magnetic flux </>) we get v a in terms of Vd as 
follows: 

v a = Tv d} T = B + e^AT (l - e^BT )- l A. (5) 
Since T is a transfer matrix, it must satisfy the current conservation condition: 

T*T = 1, with T* = (6) 

and the time-inversion invariance: 

f = T if = 0, with f = o{To x . (7) 

The same equations hold for T . (^]) and (0) constrain the form of the matrices A and B. 
Eq. (^) simply tells us that A = A and B = B. Notice that 2x2 matrices such that X = X 
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have the form X = xq + x ■ a, where a is the vector of Pauli matrices, x = (xi,X2, ix?) and 
all the Xj are real. For matrices of this kind the ^-conjugation just amounts to a "space" 
reflection X* = x — x ■ a. Using these properties, and the fact that (|]) and (|7]) hold for T , 
we get the following inversion formula: 

(1 - e^BTo)- 1 = — (1 - e^T*B*), (8) 

where 

A = det(l - e^BT ) = 1 - Tr (BT ) + e 2i * det B; (9) 
thus (P) can be rewritten as 

T = -£ (AB + e^AToA - e 2i,j> AB*A) . (10) 

We now plug (|T0| ) into the unitarity condition (^), which must be satisfied by arbitrary 
values of <fi and To, and is thus equivalent to the vanishing of a matrix- valued trigonometric 
polynomial. Vanishing of the sin 2(/>-term gives 

A*B = —B*A; (11) 

for = and from the arbitrariness of To we get [] 

det A = 1 - det B. (12) 

The remaining terms give no additional conditions. 

We remark that (g) and (|7|) are also satisfied if we impose on A and B the equations which 
are obtained from (0) and (|12D by reversing the signs of the right hand sides. However, 
only the positive choice of the sign respect invariance with regard to the interchange of A 
and B. This invariance has to be satisfied for the following reason: if we connect b and d 
instead of h and c we would obtain ([TT]) and (|T2j ) with A and B interchanged. It is easy 
to verify that the new equations would not be compatible with the previous ones if we had 
chosen the negative signs. 

Using (|TTD and ([12]) we can rewrite T in the simpler form 

T= ^-[(2cos0- Tr£T ) ■ B + AT A\. (13) 



One can easily verify that the general solution of (|TT| ) and ( |I2"D depends upon six real pa- 
rameters and is given by 

B = kSA (14) 

with 

S = i{*- J), ±l=g 2 -M 2 , P = ±(l/detA-l), (15) 



5 Eq. 12 can also be rewritten as B*B + A* A = 1. 



3 



where A = a® + a ■ a is an arbitrary matrix satisfying A = A and the positive signs in ([T 
have to be chosen when < det A < 1 . 

Landauer's formula ([I]), using (|13D, now reads 



G -2^|QP-|A| 2 ' (16) 

with 

Q = (2cos0 - Tr BT ) ■ B 22 + {AT A) 22 . (17) 

3 Appearance of the AAS effect 

We shall consider the geometry depicted in Fig. 1, where the loop of length L is constituted 
by a 1-dimensional disordered conductor. In this case the transfer matrix describing the 



propagation through the disordered loop can be parametrized as follows pT 



T _ I a P \ _ ( e i{ - as+kL ^ cosh V e^^+^sinhT \ , , 

°~ \{3* a* ) ~ { e^ +fcL )sinhr e-^ +fcL Whr J ' ^ } 

where a s , (3 S and Y are functionals of the impurities potential and can be considered as slow 
variables (see |12|] for details). 



The main point in the approach of Refs. [0, |], |K], [TT|] is the averaging over an ensemble 
of impurities potentials. However, in the case of a given sample we also have to separate 
fast-phase averaging and ensemble averaging. Fast-phase averaging can be performed e.g. 
in the following way: 

1 r L+AL 

G L = — GdL' (l/ife<AL«4 (19) 



AL Jl 

where i is the localization length. This average must be considered since the length L is 
defined with a precision AL of the order of the size of the intermediate domain where the 
wave function is essentially three-dimensional. From (fl8|) it follows that only functions of 
| a | 2 and \/3\ 2 are non- vanishing after fast-phase averaging. 

We are now in a position to prove the appearance of the Al'tshuler-Aronov-Spivak effect 
in our model. From (|^) and ( |T7| ) it follows that the coefficients of cos 0-terms appearing in 
the r.h.s. of (fIB| ) are linear combinations of a, a*, (3 and (3* . Therefore after performing the 
average (|19D only even powers of cos <j) survive and G becomes a function of period n instead 
of 2ir. An explicit example is given in the following section. This shows that the conclusion 
of Ref. H that at \(3\ ^ the cos 0-oscillations of the conductance survive after fast-phase 
averaging is wrong 0. 



We must stress here that the phases of the non-diagonal elements of (|18|) are fast variables. 
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4 Explicit computation 



Let us first of all consider the simple case of the following 1-parameter solution of (|Tl|) and 

my. 



A = y/l - K 2 , B = iKCT 3 . (20) 

From (p~6|) we get 

e 2 4k 2 Im 2 a + 4k(1 + k 2 ) cos01ma + 1 + 2k 2 cos 20 + k 4 

2^fi (1-k 2 ) 2 \(3\ 2 ' (21) 



After the fast-phase averaging we obtain 

e 2 1 + An 2 cos 2 + k 4 1 
2^h (1 - k 2 ) 2 p 



_ c ± -r -±n, l-uo y -r ^ / 99 n 



G ~ " 2 2 ^ (23) 
Goo ~ 2nh(l-K 2 ) 2 - {Z6) 

In the general case the expressions of Gl in terms of |a| 2 , \f3\ 2 and cos 2 is rather cumber- 
some. If L/£ 3> 1 we have 

|/3| 2 = exp log |/?| 2 ~ exp (log |/?| 2 ) (24) 



since the logarithm of \a\ 2 and \(3\ 2 in the limit L/l ^> 1 is additive JT3|, |14[. For large L/£ 
the behaviour of (log |/3| 2 ) is well known: 

log|/3| 2 > ~ 2^. (25) 



Thus, the values of \a\ 2 and \/3\ 2 are large enough and Gl can be represented as a series 
expansion in |/3|~ 2 . At the first order we get Gl = Goo + C + D cos 2 where C and .D are of 
the order exp(— 2L/£). This is in accordance with estimations |Tj| made in the framework 
of weak-localization theory. 



5 Chain of loops 

In this section we want to show that the interference effects due to the presence of loops can 
in principle lead to the derealization of the wave function. 

The transfer matrix of a conducting loop with coefficient of internal scattering of order 
5x (i.e. with A = 1 + A' 5x, B = B' 8x) is obtained by linearizing ([13]): 

R=l+iQ5x (26) 

where G is a matrix verifying 0* = — 0, = 0. Let us consider a series of N such loops 
connected by short wires of disordered conductor of length Sx, and let us denote by Xi the 
position of the z-th loop. The transfer matrix corresponding to the z-th loop is 

Ti = l + (i (p(xi) s 3 + (+{xi) s- + (-{xi) s + )5x = 1 + Qi Sx, 
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where iip and (± are local random fields, which in terms of the random potential V(x) have 
the form (p(x) = —V(x)/k, (±(x) = ±iV(x) e ±2tkx /2k (see [12] for details). Introducing 



the unitary matrices Ui = (R) 1 we can write the transfer matrix for the loop chain in the 
following form: 

N 

T = R N HU*TiUi 

i+1 
N 

= R N Y[(l + U:Q l U i 5x) (27) 

i=l 



exp(iOL) -Vexp J d x U*(x) Q(x) U(x) 



where V denotes the operator ordering of the exponential along the integration line. Let us 
impose the resonant condition that O be diagonal. Choose k so that 

iQ = -ika 3 , (28) 

thus giving 

U(x) = e- tkxa3 (29) 

and 

U*(x)Q(x)U(x) = ^-( 1 1 (30) 

The path-ordered exponential is easily computed since the matrix ([30]) is factorized: 

T = e -i k L«J^^ k S- L V{x)dx -±fL L V(x)dx \ 

{ ±£ L V{x)dx l-&£ L V{x)dx) 1 ] 

It is worth noting that the non-exponential dependence of the matrix T on f_ L V(t) dt follows 
from the nilpotency of the matrix (p0[). The resistance of the loop chain is then 



Ti 2 T 21 = f dxdx'V(x)V{x'). (32) 

This result is exact and does not depend on taking any average (it refers of course to a 
very particular case). For a random function V(x) with space-homogeneous statistics |13 
the right-hand side of (|32|) is a self- averaging quantity. Therefore the bulk resistivity of the 
chain is equal to l/£, where I = 4k 2 /D is the localization length in the absence of loops. 
We remember that in the absence of loops the resistance growths exponentially with L. The 
linear behaviour of the resistance (Ohm's law) obtained in place of the exponential one shows 
that the localization length diverges. 
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